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We develop a linked cluster method to calculate the spectral weights of many-particle excitations 
at zero temperature. The dynamical structure factor is expressed as a sum of 'exclusive' structure 
factors, each representing contributions from a given set of excited states. A linked cluster technique 
to obtain high order series expansions for these quantities is discussed. We apply these methods to 
the alternating Heisenberg chain around the dimerized limit (A = 0), where complete wavevector and 
frequency dependent spectral weights for one and two-particle excitations (continuum and bound- 
states) are obtained. For small to moderate values of the inter-dimer coupling parameter A, these 
lead to extremely accurate calculations of the dynamical structure factors. We also examine the 
variation of the relative spectral weights of one and two-particle states with bond alternation all 
the way up to the limit of the uniform chain (A = 1). In agreement with Schmidt and Uhrig, we 
find that the spectral weight is dominated by 2-triplet states even at A = 1, which implies that a 
description in terms of triplet-pair excitations remains a good quantitative description of the system 
even for the uniform chain. 

PACS numbers: 75.40.Gb, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 



Controlled and systematic calculation of dynamical properties of quantum lattice models remains a challenging 
computational task. Despite much recent progress in developing computational methods, such as the density matrix 
renormalization group, quantum Monte Carlo and series expansions, the dynamical properties, especially those as- 
sociated with multiparticle excitations, remain poorly understood. In many systems, these multiparticle excitations 
have relatively small spectral weights. But in low-dimensional systems, they can become extremely important and 
even dominate the spectral functions. The increased frequency and wavevector resolution of various dynamical experi- 
ments, especially neutron scattering, necessitates going beyond the single-particle picture and obtaining quantitatively 
accurate results for the full dynamical structure factors. 

One effective way to study quantum lattice models has been by using high order power series expansions in a suitable 
coupling constant. These methods have long been used to study ground state properties and elementary excitation 
spectra. Recently, they have been extended to multiparticle excitation spectra as well [1-3]. In particular, the linked 
cluster expansion method [4] is a particularly efficient way to carry out these calculations, which provides substantial 
internal checks on the validity of the expansions, and allows one to carry out the calculations by automated computer 
programs. 

Here, we develop a general linked-cluster formalism to calculate the single-particle and multi-particle contributions 
to the dynamical structure factor. We apply the method to the alternating Heisenberg chain (AHC), where expansions 
are done around the strong coupling limit of decoupled spin dimers up to 14th order in the ratio of coupling constants. 
We calculate various properties of the multiparticle continuum and bound states for a number of parameters including 
those appropriate for the material Cu(./V 03)2 -2. 5D20, which has recently been studied by neutron scattering [5,6]. 

The alternating Heisenberg chain is an excellent test-bed for our present purposes. It is a simple isotropic spin 
system with a finite energy gap (see Section III). The ground state is a spin singlet, and the lowest excitations are S 
= 1 triplet states. Above this band of 1-particle states there lies a 2-particle continuum. It was noted recently [7-9] 
that there are also 2-particle bound states with both 5 = and S = 1 below the continuum, as well as quintuplet 
antibound states above it. Trebst et al. [1,10] found, in fact, that over a wide range of parameters there exist two 
singlet and two triplet bound states near the Brillouin zone boundary kd = n (where d is the inter-dimer spacing). 
Hence the model displays some interesting multi-particle dynamics which can be explored both theoretically and 
experimentally. Neutron scattering experiments, however, will only be sensitive to the triplet bound states, which is 
what we focus on here. 

There has been much discussion in the literature [11-13] about the behaviour of the alternating chain as it approaches 
the uniform limit, which is a critical point of the model. The uniform chain is gapless, and is known to exhibit 
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"spinon" excitations. The alternating chain is gapped, and is described in term of triplet excitations. The debate has 
concerned the behaviour of the model near the uniform limit, and how the transition is made from the "triplet" mode 
of description to the "spinon" mode. We discuss these questions further in Section IV. 

We also directly obtain series expansions for sum- rules representing the total contributions of two-particle excitations 
to the dynamical structure factor summed over all wavevectors. Comparison of these with the static structure factor 
and the one-particle spectral weights, and extrapolations using approximant methods, leads to the conclusion that 
just keeping the one and two-particle excitations leads to a highly accurate description of the full dynamical structure 
factor, and this description remains quantitatively valid even in the uniform chain limit. 

The plan of this paper is as follows. First, we discuss the formalism for calculating spectral weights using the linked 
cluster method. This is followed by detailed calculations of the spectral functions for the alternating Heiscnberg 
chains. These arc followed by our conclusions. 

II. FORMALISM 

We follow the formalism of Tennant et al. [6] . The inelastic neutron scattering cross-section [14] 

d 2 a ..x-^ kf , 
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is proportional to the "dynamical response" S a/3 (k,uj) where k is the wavevector transfer, F(k) is the magnetic form 
factor, N is the number of scattering centres, fcj and kf are the momenta of initial and final neutron states respectively, 
and a = x,y,z are Cartesian spin coordinates. The dynamical response is the space and time Fourier transform of 
the spin-spin correlation function 

S^(k )W ) = — 2 J jxp^t + k-in-r.MSfi^mdt (2) 

where i and j label sites of the system. At temperature T = 0, the dynamical response becomes 

1 f°° 

S^(k,uj) = jM^t + k- (r, - rj ))]^o|5;(t)Sf(0)|^)dt (3) 

where |^o) is the ground state of the Hamiltonian. This quantity is referred to as the neutron scattering "structure 
factor" . For the alternating Heiscnberg chain, spin conservation and isotropy in spin space ensure that S a ^ = for 
and all diagonal spin components are equivalent, 

S xx (k, u) = S vv (k, uj) = S zz {k, lo), (4) 

and 

S xx {k,u)= l -S- + {k,w). (5) 
Henceforth we concentrate our attention on S~ + (k,ui). 

A. Integrated Structure Factor 

Integrating equation (2) over energy, we get: 

/oo 
duj,S- + (k,L0) 
-oo 

= ^Eexp[zk.(r i -r j )](575+) (6) 
which is just the Fourier transform of the spin-spin correlation function at t — 0. By translation invariance, 
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(SrS+) = C- + (i*,r i -r j ) (7) 

where i* labels the position of i within the unit cell. In the present case, the unit cell consists of a single dimer. 
Therefore 

5 _+ (k) =^C- + (i*,£)cxp[ik-<5] (8) 
s 

This quantity we shall refer to as the integrated structure factor. 

B. Exclusive Structure Factors 

Inserting a complete set of eigenstates of H between the spin operators in equation (3), we can express the 
spin structure factor as a sum over "exclusive" structure factors 

S-+(k,c)=]T^ + (k,c) (9) 

A 

where 

sx + (Kuj) = —Y, J JteM^t + K-^-T.MMSjm^^^s+m^) (io) 

Each exclusive structure factor 5^ + (k, u>) gives the intensity of scattering from \ip ) to a specific triplet excited state 
IV'a) [15]- In the Heisenberg picture S~(t) = cxp(iHt)S~ (0) exp(— iHt) gives trivial exponentials in t, so the time 
integral simply gives an energy-conserving delta function 

SX + (K^) = ^<5( W -S A + So)|^(^ A | > S+|Vo)exp[zk.r J ]| 2 (11) 

i 

Assuming the states \tpA) are eigenstates of momentum, the matrix elements of the spin operators at translationally 
equivalent sites are equal modulo a plane wave 

(V>A(k)|S+|Vo> = (VA(k)|S+|Vo)expHk • (r, - i>)] (12) 
Then the sum over all sites i can be reduced to a sum over sites i* in the unit cell 



N 2 i 2 

5^+(k, w) = -j£-5(u -E A + E )\ ^ A (k)|S+|^ ) exp[zk • r>] (13) 

i* 



where N c is the number of unit cells on the lattice. It is convenient to define the "reduced exclusive structure factor" 
as 

SX + (Kuj) = N c S(co -E A + Eo)\ ^(VA(k)|5+|^ }exp[ik • r 4 ,]| 2 (14) 

i* 

If we sum over momenta, we obtain the autocorrelation function 



k 

made up of exclusive contributions 

*aM4E S a + (M) (16) 



27r , 



For the 2-particle continuum, one can also obtain the total auto-correlation function $a = ^ w $a (<*>). 

We turn now to a discussion of algorithms for the calculation of exclusive structure factors within perturbation 
theory. Efficient linked cluster expansion methods have long been known [16-18,4] for calculating bulk properties of 
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a quantum lattice system. Similar methods for the calculation of 1-particle spectra were developed by Gelfand [19], 
and were extended to 2-particle spectra by Trebst et al [1,2]. 

Series methods for the calculations of exclusive 1-particle structure factors or spectral weights have been developed 
previously [20], and have been applied in several places (e.g. [21,9,22]). In this formalism, a cluster expansion was 
carried out directly for the structure factor itself. This formalism is inapplicable to the 2-particle bound states, 
however, because one needs to know the wavefunction of these states beforehand. Some leading order hand calculations 
for 2-particle states have recently been made [15,6]. The method proposed here performs a cluster expansion for 
the 'exclusive matrix elements', solving the effective 2-particle Hamiltonian to generate the required wavefunctions 
implicitly. A short paper giving some of our results has appeared recently [23] . We also note at this point that Knetter 
et al. [3] have used an alternative approach based on 'continuous unitary transformations' which is also capable of 
giving bound state energy spectra and structure factors to high order and in great detail, but this approach is more 
suitable to low dimension case. 

Let us suppose that the Hamiltonian can be decomposed 

H = H + W (17) 

where Hq is the unperturbed Hamiltonian and V is to be treated as a perturbation. We aim to expand the multiparticlc 
dispersion relations and structure factors in powers of the parameter A. For illustrative purposes, we shall use the 
language of the AHC, but the formalism can be applied more generally. 



C. 1-particle states 

At zcroth order, the 'single-particle' excitations in this model consist of triplet excitations on a single dimer which 
can be labelled \i[)\(m)), where to is the position of the excited dimer, and A labels the angular momentum eigenstate 
(i.e. S 2 , S z , in this case for AHC, the label A must correspond to the (S, S z ) = (1, 1) state). In momentum space, the 
eigenstates will be 

|Va(£0) - -= ^exp(ik • r m )|^ A (m)) (18) 

» c m 

This labelling can be retained at higher orders in perturbation theory as A is raised from zero. 
Then the matrix element 

(V> A (fc)|^|Vo> = -^^expHkT ro )<V>A(m)|S+|V>o> (19) 

* c m 

It follows from translation invariance that the matrix element {i()\(m)\S* \i/)o) in this expression is a function of 
(i-j — r TO ) only, i.e., we can define the exclusive matrix elements 

n^{i*,S)=(Mm)\S+\^ ) (20) 

where 6 is the distance between i and to, S = — r m , and i* labels site i within the unit cell as before. 
The reduced exclusive structure factor is 

^ + (k) = \ Y / nA(t*,S)exp[ik-8}\ 2 (21) 

i*,6 



D. 2-particle states 



Since two triplets can combine to give total spin 1, there will also be non-zero spin structure factors for 2-particlc 
states. The 2-particle states can be labelled according to their unperturbed counterparts in position space \i/)\(m, n)) 
where to, n label the two dimer positions, and A the corresponding total 2-particle angular momentum states. Eigen- 
states of the Hamiltonian can then be expanded 

|^A(k)) = ^ E A(r m r n y k <^+^\yj A (m, n)) (22) 

v c m,n 
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where f A (r m — r„) is the 2-particle wavefunction, which only depends on the relative distance r m — r„. 
Similarly one can define the 2-particle exclusive matrix elements 

n 2 A p (i;m,n) = W A (m,n)\S+\ik) (23) 

The translation invariance implies that Cl A F '(i;m,n) only depends on the relative distance between to, n and i, i.e. 

Q A p (i;m,n) = Q A p (i*,r,S) (24) 

where r = (2r; - r m - r„)/2 and 8 = r m - r„. 

Then, we can get the reduced 2-particle exclusive structure factor as 

5A + (k,c^) =6(lj-E a + Eo)\ ]T tf A p (i*,r,S)f A (6)exp[ik-r}\ 2 (25) 

i* ,r.d 



E. Cluster expansion 

It is easy to see that the exclusive matrix elements il A p (S) and A p (r, S) obey a simple 'cluster addition' property. 
If a cluster C is made up of two disconnected sub-clusters A and B, then 1 

(26) 

fi A p ' C (r, S)= n 2 A p - A (r, 6) + n 2 A ^ B (r, 6) (27) 

where f2 A p ' A (£) is trivially zero if cluster A does not contain both dimer to and site i, by conservation of spin, and 
similarly for 0^ p ' A (r, 5). It follows that the elements Q A admit a linked cluster expansion 



where the sum over 7 denotes a sum over all connected clusters. Correspondingly, the perturbation series expansion 
for fl A could be formulated in terms of a diagrammatic expansion where only connected diagrams contribute, although 
we will not elaborate on this approach here. 

An efficient linked cluster algorithm for calculation of the structure factors can now be formulated, following Trebst 
et al. [1,2]: 

i) Generate a list of connected clusters 7 appropriate to the problem at hand (in the present case, they will simply 
consist of chains of dimers of different lengths); 

ii) For each cluster 7, construct matrices for the Hamiltonian H and spin operators Sf in the basis of singlet and 
triplet dimer states corresponding to H n ; 

iii) 'Block diagonalize' the Hamiltonian by an orthogonal transformation 

H ctf = O t HO (29) 

as outlined by Trebst et al [1,2], constructed order-by-order in perturbation theory so that the 1-particle states 
sit in a block by themselves, and similarly the 2-particlc states, etc; with this one can compute the exclusive 
matrix elements ft A . 

iv) Subtract all sub-cluster contributions to get the cumulant f2A(£); 



1 We suppress the starred cell index i* , henceforth. 
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v) Insert the cumulant Oa(5) in Eq. (28), and hence one can build up the exclusive matrix elements for the bulk 
system. 

vi) For the 1-particle case, one can insert the exclusive matrix elements into Eq. (21) to get the series for the 
exclusive structure factors. For the 2-particle case, one still needs to solve the effective 2-particle Hamiltonian 
to get the wavefunctions /a for the possible bound states and continuum: this can be done by using the finite 
lattice approach [2] . To get the series solution rather than the numerical solution for the wavefunctions /a and 
the 2-particle dispersion relation, one has to use degenerate perturbation theory. With the results for exclusive 
matrix elements f^A and wavefunctions /a in hand, one can get the exclusive structure factor S\ through Eq. 
(25). 

III. RESULTS FOR THE ALTERNATING HEISENBERG CHAIN 

We apply this method to investigate the spectral weights of the alternating Heisenberg chain, which can be described 
by the following Hamiltonian 

H = ^ S 2i • S 2 i+i + AS 2 i-i • S 2 i (30) 

where the Si are spin-i operators at site i, and A is the alternating coupling. Here we assume that the distance 
between neighboring spins are all equal and the distance between two successive dimers is d, and assume d = 1 if it 
does not appear explicitly (note this is different from our previous paper [10], where we had taken the lattice spacing 
a to be 1). 

There is a considerable literature on this model, which has been reviewed recently by Barnes et al. [15]. At A = 0, 
the system consists of a chain of decoupled dimers, and in the ground state each dimer is in a singlet state. Excited 
states are made up from the three triplet excited states on each dimer, with a finite energy gap between the singlet 
ground state and the triplet excited states. This scenario is believed [24-26] to hold right up to the uniform limit 
A = 1, which corresponds to a critical point. At A = 1, we regain the uniform Heisenberg chain, which is gapless. 

Several theoretical papers [11-13] have discussed the approach to the uniform limit. Analytic studies of the critical 
behaviour near A = 1 [11] have related the alternating chain to the 4-state Potts model, and indicate that the 
ground-state energy per site e (A), and the energy gap A(A) should behave as 

e (A)-e (l)~,5 4 / 3 /|ln(<V<5o)| 

(31) 

A(A) ~ 6 2 ' 3 /y/\H6/6 )\ 

as A — > 1, where S = (1 — A)/(l + A). The logarithmic terms are due to the existence of a marginal variable in the 
model. 

Numerical studies of the model include series expansions [24,25,22], and exact diagonalizations for finite lattices 
[27,28]. Recently, Papenbrock et al. [29] have carried out density-matrix renormalization group studies on lattices up 
to 192 sites in extent. They conclude that the data for the ground-state energy and triplet energy gap are consistent 
with Eq. (31), but with surprisingly large scale factors So in the logarithms. Dlog Pade analysis of the series supports 
these conclusions [22]. 

The 2-triplet bound states were previously studied by Uhrig and Schulz [7] using an RPA approach. They found a 
singlet bound state below the 2-particle continuum for all momenta k and over the whole range of 1 > A > 0. They 
also predicted a triplet bound state and a quintuplet antibound state near kd = n for small A. These conclusions were 
supported in later studies [8,9,30]. Barnes et al. [15] and Tennant et al. [6] have shown how to calculate exclusive 
structure factors for the 2-particle states by low-order series expansions in A, and have made a comparison with 
experimental data for the copper nitrate material, Cu(iV0 3 ) 2 . 2. 5D 2 0. They find that the experimental data are 
consistent with the existence of a bound state, but do not yet constitute definitive proof of it. They highlight the 
need for more powerful perturbative techniques to calculate these multiparticle cross sections: the present paper is 
aimed at meeting this need [23]. 

We have recently made an extensive study of the two triplet bound states of this model using high-order series 
expansions [1,10]. We found that in fact there are two singlet (Si and S 2 ) and two triplet (Ti and T 2 ) bound states 
below the two-particle continuum, and two quintet antibound states (Qi and Q 2 ) above the continuum, for A not too 
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large. Meanwhile, Schmidt and Uhrig [31] have used a different technique, the 'continuous unitary transformations' 
(CUTS) method, to study the model at high orders in perturbation theory. We shall compare our results with theirs 
in what follows. 

Here we study the structure factors using series expansions. We have computed, up to order A 13 , the series for the 
integrated structure factor S(k) and the exclusive structure factors for the 1-particle triplet state Si p (k). For the 
2-particle states, we also computed the total 2-particle structure factor S 2 p (summed over all 2-particle states) up to 
order A 13 , and to order A 12 for the exclusive matrix elements f2^ p (r, <5). With these we can compute the structure 
factors for the 2-particle triplet bound states (i.e. St x for T\ and St 2 for T 2 ,), and for the 2-particle continuum 
S2pc up to order A 14 . Integrating the structure factor over momentum k, we can also compute the series for the 
auto-correlation function defined in Eq. (15). We have computed the auto-correlation functions for the 1-particlc 
(<&ip), and 2-particle sectors ( < &2 P ), and the individual auto-correlation functions for 2-particle bound states T\ (3>Ti) 
and T 2 ($t 2 ) an d the 2-particle continuum (<&2pc)- Full scries for S(k) and S\ p (k) are given in Tables I and II. Some 
other selective series (at kd — > 0, kd = it, and 2tt) are given in Table III, other series are available on request. 

The 1-particle structure factor has been computed to order A by Barnes et al. [15], but our series disagree with 
their results from second order. In their calculation, they neglect the first term in their Eq. (56); they claim this 
term is zero due to a symmetry relation, but actually this only holds in leading order. Recently Muller and Mikeska 
[32] have extended the series for the 1-particle structure factor S\ p (k) to order A 10 , and our series agree with theirs. 
The auto-correlation functions for 1-particle ($i p ) and 2-particle ($2 P ) states have been recently computed up to 
order A 7 by Schmidt and Uhrig [31], and our results also agree with theirs, but extend the series by up to 6 terms. 
As a byproduct of our calculation, we have computed the series for 1-triplet excitation spectra up to order A 13 , and 
the series for 2-triplet excitation spectra up to order A 12 , this extends the previous calculations [10] by two terms for 
1-triplet excitation spectra, and by 1 term for 2-triplet excitation spectra. These series are available on request. We 
use Dlog Pade approximants and integrated differential approximants [33] to obtain numerical results up to A = 1. 

A. The integrated structure factor S 

The integrated structure factor 2 S has been computed up to order A 13 . It can be expressed as 

oo 2m+l 

S(kd) = 1 - cos(fcd/2) + a n,mX m cos(nkd/2) (32) 

m— 1 n— 1 

The series coefficients a n , m are given in Table I. Note that there are two general relations for these coefficients: 

2m+l 
n=l 

(33) 

With this, one can easily prove that S(k) can be written in the following form with a common factor sin 2 (kd/ 4) at 
all orders 

oo 2m 

S{kd) = sin 2 (jfed/4) Yl VmA m cos(nfcd/2) (34) 

where one of the coefficients 62m-1.n1 happens to be zero always for any m, due to the second relation in Eq. (33). 
For example, the series up to order A 4 can be written as 

kd [ A 2 / kd 

S{kd) = sin 2 (— ) 2 + A cos(kd) + — 2 + 4 cos( — ) + 3 cos(fcrf) + 3 cos(2 kd) 



4 
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2 We drop the superscripts — h henceforth. 
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A 3 / kd 3 kd 

+— 12 + 24 cos( — ) + 13 cos(fcrf) + 16 cos( — -) + 28 cos(2 kd) + 15 cos(3 kd) 
96 y 2 2 

A 4 / kd 3kd 5 kd 

+— — 346 + 692 cos(— ) + 669 cos(kd) + 912 cos(-— ) + 784 cos(2 kd) + 304 cos( — -) 
4o(J8 V 2 2 2 

+852 cos(3 kd) + 315 cos(4 kd)) + 0(A 5 ) (35) 

To analyze the series S(kd), we firstly consider S(kd) near two special momentum points, kd = and kd = 27r. We 
know that S = at k = 0, but it is interesting to study how it vanishes as k — > 0. For small A, it is trivial to see from 
Eq.(35) that S oc k 2 as k — > 0. Thus we can define i? = linifc^o S '(kd) / '(kd) 2 , and the series for R is given in Table 
III. Applying Dlog Pade approximants to this series, we find that R diverges at A = 1, with a critical index about 
-0.71 3 . This implies that for the uniform chain, A = 1, S(kd) no longer vanishes quadratically with k. Naively one 
might then expect it to vanish linearly as k — > 0: some evidence for this is shown in Fig. 1. 

For the uniform chain A = 1, and near kd — > 2n 7 Affleck has argued [34] the asymptotic form for S(kd) as 

^ (fcd) = 3727)^ |ln(7r_M/2)|3/2 (36) 

This implies that for kd = 2n and as A — > 1, the asymptotic form for S(2n) diverges as 

S(2ir) oc [- ln(l - A)] 3 / 2 A -> 1 (37) 

If so, one should be able to see a critical point at A = 1 with critical index -1 if one applies Dlog Pade approximants 
[33] to the series dS 2 / z /d\. This is indeed the case as we can see from Table IV. Note, however, that very similar 
results are obtained for any index between 3/2 and 1, so this is not a very sensitive test. 

Assuming the above asymptotic form for S(2tt) is correct, we can also estimate the prcfactor of this asymptotic 
form by using integrated different approximants [33] to extrapolate the series for S(2n)/[1 — ln(l — A)] 3 / 2 to A = 1. 
The result is: 

S(2i:) = 0.19(2)[- ln(l - A)] 3 / 2 A -» 1 (38) 

The prefactor agrees with Affleck's prediction of 8/3(2tt) 3 / 2 = 0.16932. 

For < kd < 2tt, one expects S to be finite for any A. The results for S versus momentum k for A = 0, 0.6, and 1 
are shown in Fig. 1. Note that f Q n S(k)dk — 2ir (here we set d = 1), independent of A, so the area under each curve 
is the same. Also shown in the figure are the results for 5" = 6S[— 27rln(l — ^)/k]~ 3 / 2 at A = 1. 

For fixed values of k, Fig. 2 shows the integrated structure factor S versus A, where for each value of k, about 20 
different integrated differential approximants to the series are shown. We can see that the results converge very well 
out to A = 1. The logarithmic divergence as A — > 1 for the case kd = 2n is clearly evident. 

B. The 1-particle spectral weight 

The exclusive 1-particle structure factor Si p (kd) has been computed up to order A 13 . 5i p can also be expressed as 

S lp (kd) = a«, m A m cos(nkd/2) (39) 

The series coefficients a n , m are given in Table II. The coefficients o„ jm also satisfy Eq. (33), and Si p (kd) also has a 
common factor sin 2 (fcc?/4) for all orders. The series up to order A 4 can be written as 



kd 

Sip(kd) = sin 2 ( — ) 2 + A cos(fcrf) + ^- ( -1 + 8 cos(^) + 2 cos(kd) + 3 cos(2 kd) 



A 2 / . . „ ,kd, 
~2 



A 3 / kd 3 kd 

— 36 cos( — ) + 17 cos(fcrf) + 8 cos(— ) + 32 cos(2 kd) + 15 cos(3 kd) 



3 Within the accuracy of our calculations this is consistent with a value of 2/3 with logarithmic corrections 
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0.5 1 1.5 2 

kd/7T 

FIG. 1. The integrated structure factor S versus fe for A = (dotted line), 0.6 (dashed line), 1 (crosses). Also shown is the 
quantity S' = 6S[-2n ln(l - A)/fe]- ; V 2 for A = 1 (squares). 




FIG. 2. The integrated structure factor 5 versus A for kd = tv/2, it, 3tt/2 and 2-k. For each value of fed, about 20 different 
integrated differential approximants to the high-temperature series are shown, though most of them are indistinguishable on 
the scale of this figure. 
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A 

FIG. 3. The relative 1-particle weight Si p /S versus A for kd = 0+, 7r/2, 7T, 37r/2 and 2-7T. Several different integrated 
differential approximants to the series are shown. 



A 4 / kd 3 kd 5 kd 

+ — — ( 342 + 384 cos( — ) + 784 cos(fcd) + 856 cos( — -) + 675 cos(2 kd) + 368 cos(— -) 
4608 V 2 2 2 

+804 cos(3 kd) + 315 cos(4 kd)) + C>(A 5 ) (40) 

Let us analyze the behaviour of Sip when A — * 1. For kd = 2ir, our analysis shows that Si p is finite as A — + 1, so 
Sip/S vanishes like l/[-ln(l - A)] 3 / 2 . For the case kd = 0, we again define Ri p = lim^^o S\ p (kd) / (kd) 2 (note that 
Rip differs from R from order A 4 ). The Dlog Pade approximants to Ri p show that it diverges at A = 1 with critical 
index about -0.71. This again implies that for the uniform chain case (A = 1), Si p no longer vanishes as k 2 . Now let us 
analyze the series for the relative spectral weight R lp /R. The Dlog Pade approximants to the series Ri p /R show no 
singularity near A = 1. This means that as A — > 1, i? lp /i? remains finite. Our extrapolations show i? lp /i? = 0.993(1) 
at A = 1. 

For < kd < 2ir, the analysis of the series Si p by the Dlog Pade approximants shows that it vanishes with a 
behavior close to (1 — A) 1 / 3 . Since S remains finite, we thus expect that Si p /S vanishes like (1 — A) 1 / 3 . This agrees 
with the analysis of Schmidt and Uhrig [31], who argued that the 1-particle spectral weight should vanish like y/A, 
i.e. like S 1 / 3 /] ln((5/ <5o) | x ^ 4 , where 8 = (1 — A)/(l + A). Fig. 3 shows the relative 1-particle weight S\ p / S versus A at 
selected values of kd. It can be seen that for any non-zero value of k, Si p /S decreases abruptly to zero as A — > 1. 
Only at kd = 0+, does Si p /S remain finite (about 0.993) in the limit A = 1; but by then S has itself decreased to 
zero. 

The overall picture that emerges is that as A — * 1 the triplet energy gap goes to zero, and the spectral weight 
associated with the 1-triplet state also vanishes [22]. This would seem to agree with the idea that the "spinons", 
rather than the triplet states, act as the elementary excitations in the model in the uniform limit A —* 1. 

Fig. 4 shows the relative 1-particle weight S lp /S versus k for A = 0.27, 0.5, 0.7, 0.8, 0.90367, where A = 0.90367 
is the estimated coupling for CuGeOa [7]. For A = 0.90367, we can see that Si p /S is about 1 for small k, and has a 
minimum around kd = 1.6n with value 0.43, whereas Uhrig and Schulz [7] estimated Si p /S = 0.28. 



C. The 2-particle spectral weight 



As we mentioned before, there are two triplet bound states (T\ and T2) in this system. Before we discuss the 
spectral weight for individual states, let us first discuss the total spectral weight (the sum rule) for all 2-particle 
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FIG. 4. The relative 1-particle weight S lp /S versus k for A = 0.27, 0.5, 0.7, 0.8, 0.90367. The results of the three highest 
orders are plotted. 



states. 



1. The total 2-particle spectral weight 

The total 2-particlc structure factor (S^p) has been computed up to order A 12 . This quantity S 2p (kd) also has a 
common factor sin 2 (fc<i/4) for all orders, and the series up to order A 4 can be written as 

^ /, ,n . ?,W,[,9 . A,kd. A 3 / „ ,3kd.\ . o,kd. A 4 / „ .kd. „ ,„ ., ,. 

S 2p (kd) = sin 2 ( T ) [A 2 sin 4 ( T ) + — (3 + 2 cos( — )J sin 2 ( T ) + — ( - 57 + 372 cos( — ) - 143 cos(fcd) 

3 fed 5 kd \ 1 

+64 cos( — ) + 108 cos(2 kd) - 64 cos( — ) + 48 cos(3 kd)) + 0(A 5 )J (41) 

Here S^p oc k 2 in the limit k — > 0, and so we again define i?2 P = limfc_o S 2p (kd) / (kd) 2 . Note that i?2 P is nonzero 
from order A , while S 2p with k ^ is nonzero from order A 2 . 

Fig. 5 shows the total 2-particle weight S 2p /S versus k for various A. For large A, S 2p /S has a maximum around 
kd = 1.6tt. For A = 0.6, it has a value of about 0.1848 at its maximum, while the value for k — is only 0.00177, 
about 100 times smaller than the maximum. 

For fixed values of k, Fig. 6 shows the relative total 2-particle weight S^p/S* versus A for kd = ir/2, ir, 3n/2 and 27r. 
We can see that there are some sharp increases near A = 1 which make it difficult to estimate the results; nevertheless, 
we estimate, at A = 1, the 2-particle states have about 90% of the total weight at kd = 3n/2, about 80% of the weight 
at kd = 2ir, about 65% of the weight at kd — n, and only about 10% of the weight at kd — n/2. We also show the 
results for 155 , 2 P /(A 2 5') in the special case k = 0, which show quite different behaviour from other k in this figure: 
as A increases, it increases firstly, then decreases, so that the relative weight at A = 1 is tiny (less than 0.5%). This 
agrees with the fact that the relative weight for the 1-particle state, Ri p /R, is almost 99.27% at A = 1 for this value 
of k. 

Subtracting the weight of 1 and 2-particle states from the total weight, one can obtain the remaining weight 
(S — Si p — S2 P )/S for states of more than two-particles. The results for various A are shown in Fig. 7. For small A, 
the remaining weight increases as k increases, while for larger A, it develops a peak near kd = 1.5n, and the remaining 
weight near kd = 2tt decreases as A increases. These weights are tiny, as can be seen, but increasing with A. 
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kd/n 

FIG. 5. The relative 2-particle weight S^p/S 1 versus k for A = 0.27, 0.4, 0.5, 0.6. The results of the three highest orders are 
plotted. 




12 



0.5 1 1.5 2 

kd/7T 

FIG. 7. The remaining spectral weight for states of more than 2-particles ((5 — 5i p — S2 P )/S) for A = 0.15, 0.27, 0.35, 0.4, 
0.45. The results of the three highest orders are plotted. 



2. The individual weights for triplet bound states and continuum 

For the bound state T\ , we can obtain an analytic expression for its structure factor up to order A 3 

kd A 3 / kd 
S Tl {kd) = -A 2 (1 + 2 cos(fcd)) sin 6 (— ) - — ( - 60 + 101 cos{ — ) - 102 cos(fcd) 

1H 5 krl \ kd 

+69 cos(^-) - 48 cos(2 fed) + 22 cos(^-)J sin 4 ( — ) + 0(A 4 ) (42) 

where the A 2 term agrees with that obtained by Barnes et al. [6] 

As a byproduct of our calculations, we also get an analytic expression for the coherence length [2] which is defined 

as 

T Irlf 2 

L = (43) 

where / r is the amplitude (the eigenvector) for two single-particle excitations separated by distance r (see Eq. 22). 
The result up to order A 2 for bound state T\ is 



j \ )2 

- — = -1 - 2 cos(fcd) + - [-24 - 36 cos(kd) - 17 cos(2 kd) - 5 cos(3 kd)] + — 
Ltx 2 32 



-3410 

5937 cos(fcrf) - 3768 cos(2 fed) - 1619 cos(3 fed) - 482 cos(4 kd) - 100 cos(5 kd)] + 0(A 3 ) (44) 



With these expressions, one can determine the critical momentum k c where S+1 or 1/Z/Ti vanishes. For both 5t x 
or 1/ , one gets the same k c as 

f 2tt/3 + 5A/(4V3) - 757A 2 /(192V3) + (9(A 3 ), k c d < tt . . 

c \ 4tt/3-5A/(4V3) + 757A 2 /(192V3) + 0(A 3 ), k c d > tt [ ' 

This agrees with previous results obtained from the two-particle binding energy [10]. The expressions (42) and (44) 
are valid within these regions of momentum. In the limit fc — > fc c , the behaviours of Sr^k) or 1/Lti are 
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FIG. 8. The relative weight for bound state Ti (S Tl /S) for A = 0.27, 0.4, 0.5, 0.6, 0.7. The results of the three highest orders 
are plotted. 



(k - k c )d\ 2 (24 + 29 A) /(512 V3) 
9V3{k c -k)d (24 -31 A) A 2 /512 



0(A 4 ), 
0(A 4 ), 



k c d < it 
k c d > 7r 



(46) 



l/L Tl = 



\k-k c \d (144 +12 A -443 A 2 ) 
48 V3 



(47) 



So as k — > k c , (k) and 1/Lt ± are proportional to (k — k c ), whereas the binding energy is proportional to (k—k c ) 2 
[10]. 

Integrating in Eq. (42) over the momenta given in Eq. (45), one can get the auto-correlation function for Ti as 



* Tl = A — 



1\ 23 V3 \ 3 / 433 361 V3\ 

6 + 576 ^^56^ + { ] 



(48) 



Since the bound state T2 only appears at kd = tt in the small A limit, we cannot get a similar analytic expression for 
it, nor for the 2-particle continuum. But the series for exclusive 2-particle matrix elements ftj^ have been computed 
up to order A 12 , and with this one can compute numerical results, for any given value of k, for the exclusive structure 
factor of the 2-particle triplet bound states (Ti and T2) and 2-particle continuum through the finite lattice approach 
[10]. We can also get some series in A for the structure factor over those momenta k where the bound states appear 
in the limit A — > 0, and the series for kd = tt are given in Table III. 

The relative spectral weights for Ti and T 2 versus k are shown in Figs. 8 and 9 for several values of A. The weight 
for Ti is nonzero only over a finite range of momentum where the bound state exists. In the limit A — > 0, St ± /S has 

a maximun at kd = 4 arccos (^/(5 — = 1.1177T, and as A increases, the maximum position kod moves towards 

to k n d = tt. The results for kod as function of A are given in Fig. 10. The relative spectral weight for T 2 exhibits an 
interference zero near kd ~ tt, which can be traced back to the bound-state wave function for A — > 0, where the two 
triplets are separated by an odd number of dimers [10]. The position of this zero shifts to smaller k as A increases, 
and is also shown in Fig. 10. Another interesting feature here is that the spectral weight vanishes like (k — k c ) 2 as 
k — ► k c , rather than (k — k c ) as we have seen in Eq. (46). Also shown in Fig. 10 is the location of the maximum 
relative weight for T 2 . 

The remaining weight for the 2-particle continuum S^pc is shown in Fig. 11. These curves show a dip near kd = tt, 
because the bound states absorb some of the weight in that region. 
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FIG. 9. The relative weight for bound state T2 (St 2 /S) for A = 0.27, 0.4, 0.5. The results of the three highest orders are 
plotted. 




FIG. 10. The momentum kod/n or k c d/ix where 7\ (T2) has its maximum (respectively zero and maximum) relative spectral 
weight graphed versus A. 
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FIG. 11. The relative weight for the 2-particle continuum (S2 pc /S) for A = 0.27, 0.4, 0.5, 0.6. The results of the three highest 
orders are plotted. 



D. Complete Dynamical Structure Factor 



The complete dynamical structure factor for A = 0.27 was presented in our earlier paper [23]. The dynamical 
structure factor for the 2-particle continuum at A = 0.5 is given in Fig. 12. A notable feature is the spike which 
develops as a bound state enters the continuum, discussed more fully in Ref. [23]. We will be happy to provide the 
complete dynamical structure factor for the AHC model for any value of A, if requested. 

Given the structure factor, one can also compute, for 2-particle states, the energy centroid (u>) and the width Au 
which are defined as 



J S2(u>, k)u>du> 



J 52 (w, k)duj 
Acu = (cu 2 ) - (cu) 2 



(49) 
(50) 



where 5 2 (w, k) is the structure factor of the 2-particle states, including both the 2-particle bound states and 2-particle 
continuum. These quantities can be extracted readily from experiments. Our results are given in Figures 13 and 14 
for various A. We can see that Aw has a nonzero minimum near kd = n, reflecting the presence of strong 2-particle 
bound states and a weak and narrow 2-particle continuum. 



E. Auto-correlation function 



Finally, we discuss the results for the spin auto-correlation functions. First, let us discuss the critical behaviour as 
A — ► 1. Schmidt and Uhrig [31] argued that the critical behaviour for the total auto correlation function (summed 
over u>) of the 1-particle state $i p and 2-particle states $2p should be 

$i p cx (1 - A) 1 / 3 (51) 
$ 2p cx const. + 0((1 - A) 1/3 ) (52) 

modulo logarithms, as for the structure factors. 

To verify these behaviours, we have applied Dlog Pade approximants to the series, and the results obtained from 
unbiased and biased Dlog Pade approximants are shown in Table IV. For <£>i p , the results are consistent with Eq. 
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FIG. 12. The structure factor (shifted by n/100) for the 2-particle continuum versus energy uj at A = 0.5 and kd = 27rn/20, 
n = 0, 1, 2, • ■ • , 20, Also shown as a bold solid line is the dispersion relation for the triplet bound-state T\. 
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FIG. 13. The 2-particle energy centroid (cu) versus k for A = 0.27, 0.4, 0.5, 0.6. 
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FIG. 14. The 2-particle width (uj 2 ) - {lo} 2 versus k for A = 0.27, 0.4, 0.5, 0.6. 

(51). Assuming this behaviour, we can estimate the prefactor using integrated differential approximants. The result 
is 

$i p = 1.258(2)(1 - A) 1/3 (53) 

For $ 2p , however, the results are more ambiguous. The unbiased Dlog Pade approximants to <9$2 P /<9A tend to 
exhibit defects for A < 1, and the biased critical index for d<&2 P /d\ is about -0.6, slightly smaller than -2/3. This 
might be due to an extra logarithmic correction. 

The various auto-correlation functions versus A are shown in Fig. 15, where one can see that $ lp vanishes at the 
limit A = 1, while (1 - A) _1 / 3 $i p increases almost linearly as A increases. The curve for (^lp + < f ) 2p), if we assume it 
is non-singular at A = 1 (i.e. the singularities in <j>i p and $2p cancel exactly), runs almost flat with A once we neglect 
unphysical and defective approximants: that would indicate that the 2-particle sector accounts for about 99.8% of 
the weight, even at A = 1, which agrees almost exactly with the conclusions of Schmidt and Uhrig [31]. Remarkably, 
this is much higher than the fraction of 73% for the two-spinon continuum at A = 1 calculated by Karbach et al. 
[35] from the exact solution. The result that sectors with more than 2 particles account for very little weight can 
be further understood from Fig. 7, where we can see that the remaining spectral weight near kd = 2ir for states of 
more than 2-particles actually decreases with A, for large A. Also shown in Fig. 15 is the direct extrapolation of the 
2-particle auto-correlation <J> 2p using integrated differential approximants. These extrapolations assume that there is 
no singularity in <j> 2p at A = 1, and the results give a somewhat smaller value of about 0.9 at A = 1. For the bound 
states, the auto correlation function for Ti increases, and reaches its maximum around A = 0.6, then decreases, while 
the auto correlation function for T\ continues to increase as far as we can follow it. We presume that $Ti will also 
vanish as A — ► 1. 

The autocorrelation function for the 2-particle continuum (<f>2pc(w)) is shown for various A in Figure 16. The major 
feature is a spike at w <~ 2 which becomes more prominent as A increases. This spike is due to contributions from the 
region around kd = it near the lower threshold of the dispersion curve. It becomes divergent as A — > 1, and matches 
rather neatly with a logarithmic divergence at oj = tt/2 discovered in the 2-spinon continuum by Karbach et al. [35]. 
Unfortunately, we can not compute an explicit series for <j> 2pc (o;), so no series extrapolation can be made here. 

Figure 16 also displays another small cusp or spike at lower uj, which occurs at the threshold energy where the 
triplet bound state T± merges with the continuum, i.e. where the structure factor diverges as shown in Figure 12. 
This position is marked by the arrow in Figure 16, for the case A = 0.6. As A — > 1, the T\ threshold migrates towards 
uj = 0, and the logarithmic divergence there [35] at A = 1 may be viewed as relic of the T\ bound state. 
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FIG. 16. The auto correlation lunction versus energy lu for the 2-particle continuum (<&2 pc (w)) f° r ^ ~ 0.15, 0.27, 0.4, 0.5, 
0.6. 
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FIG. 17. The average weights (Wi p and W2 P ) for 1 and 2 particle states, and ratio H / 2 P /$2p versus A. Several different 
integrated differential approximants to the series are shown. 



Note that the structure factor is dominated for A — > 1 by kd — 2tt, and so is the auto-correlation function. To get 
an idea of the multiparticle contributions to the dynamical structure factor at all wavevectors, we define the average 
relative weight, for a particular state A, as 



The series for the average relative weight for 2-particle states is given in Table III, and graphed as function of A in 
Figure 17. Once again, we see that the average weight for the 1-particle state drops to zero as A — > 1. The average 
weight for the 2-particle states has a sharp increase near A = 1 which makes it difficult to estimate the results at 
A = 1, but the result for the ratio TF^p/^p is quite flat. Hence one estimates that W2 P remains substantial, at about 
60%, as A — > 1. This is smaller than <J> 2p , showing that multiparticle excitations are more important away from the 
antifcrromagnetic wavevector kd — 2ir. 

A priori, one might have expected that the average weights for 1-particle states, 2-particle states, etc, would all 
tend to zero as A — > 1, with all the weight moving into many-particle states. Instead of that, we find that the weight 
for the 2-particle continuum remains finite and large, comparable to or by some measures even greater than that of 
the 2-spinon continuum computed by Karbach et al. [35] 



We have shown in this paper how to calculate multiparticle structure factors and spectral weights to high orders 
in perturbation theory using linked-cluster expansion techniques [1,2]. Applying these techniques to the case of the 
alternating Heisenberg chain, a detailed picture has been given of both the integrated structure factors, and the 
individual spectral weights for the 1-particle state, 2-particle bound states and 2-particlc continuum as functions of 
wavevector k. Continuing the series by means of Pade approximants or integrated differential approximants, good 
convergence is obtained from A = right up to A = 1. Hopefully, it should be possible to test these predictions against 
experiments in the near future [5,6]. 

The 1-particle energy gap and spectral weight at general momenta appear to vanish as A — > 1, following the 
behaviour predicted by Cross and Fisher [11], and already confirmed numerically by Singh and Zheng [22]. This 




(54) 



IV. CONCLUSIONS AND DISCUSSION 
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would seem to confirm the general notion that the triplets no longer form elementary excitations for the system at 
A = 1. However, the 2-triplet spectral weight remains finite in the uniform limit and, in fact, appears to form the 
major part of the total spectral weight [31]. Schmidt and Uhrig [31] already pointed out that indeed the 2-triplet states 
carry a larger portion of the total spectral weight than the 2-spinon states, calculated by Karbach et al. [35]. We also 
find that the 2-particle auto-correlation functions for triplets display singularities similar to the known singularities 
for spinons [35]. 

We find that clear precursors to many features of the two-spinon continuum for the Hciscnberg model, such as 
vanishing weight at the upper end of the continuum, divergent weight at the lower end of the continuum, cusps in 
the autocorrelation function, etc. , are already evident in the two- triplet continuum of the alternating chain. We have 
discovered an interesting result in that the logarithmic divergences at u> = in the auto correlation function of the 
uniform chain appears to be a relic of the triplet bound state T\ in the non-uniform chain. Our highly accurate 
calculations of frequency and wavevector resolved spectral-weights should prove useful in understanding the spectral 
functions of real materials. 

The crossover from elementary triplets to spin-half elementary excitations is quite different here than in the case 
of the Ji — J2 — 5 model with J2/J1 > 0.24, where the uniform limit stays spontaneously dimerized. Affleck et al 
[13] presented the following picture in the latter case. The mass gap decreases exponentially as the uniform limit 
is approached, and in that limit the low-lying spectrum consists of a gapless continuum of soliton-antisoliton states. 
Away from the uniform limit, the soliton-antisoliton states are confined by a linear potential, giving rise to a ladder 
of discrete states, which correspond to the triplet excitations and their bound states. In the triplet languange, on 
the other hand, the triplet energy gap drops to zero as the uniform limit is approached, and the multi-triplet bound 
states condense to form a continuum matching the soliton-antisoliton description. Numerical evidence confirms this 
picture [13,10]. 

The scenario appears to be quite different in the present case, where there is no next-nearest-neighbour interaction 
(J2)- Here the dimerization interaction has dimension |, giving rise to a gap which vanishes as <5 2 / 3 in the uniform 
limit. In the soliton language [36], at small but non-zero 5, the low-lying spectrum consists of a soliton, an antisoliton, 
and a soliton-antisoliton bound state, forming a degenerate triplet, plus just one other soliton-antisoliton state at \/3 
times the mass. Since the soliton and antisoliton are not confined, a soliton-antisoliton continuum occurs at higher 
energy. In the triplet language, the 1-particle triplet corresponds to the three degenerate states, and the singlet 
bound state corresponds to the higher-lying one. There is no condensation of multi-triplet bound states; instead, 
the 2-triplet continuum drops down to match the spinon-antispinon continuum. The mass and spectral weight of the 

1- triplet state vanishes as the uniform limit is approached; and naively one might expect the same to happen for the 

2- triplet, 3-triplet states, etc, leaving many-triplet states to correspond to the spinon-antispinon continuum. Instead, 
we have found that the 2-triplet states dominate the total weight in the uniform limit, in agreement with Schmidt 
and Uhrig [31]. This is a peculiar and paradoxical feature, which argues that the triplet description remains at least 
as relevant as the spinon description, even in the uniform limit. The multi-triplet states (n > 2) still appear to carry 
only a tiny fraction of the spectral weight at A = 1 . 

These findings call for further theoretical interpretation. Karbach et al. [37] have described the uniform Heisenberg 
antifcrromagnct in terms of three different elementary excitations - magnons, spinons, and 'psinons' - each of which 
can be useful in different circumstances. None of these correspond to our dimcr triplet excitations, however. Our 
work shows that a conventional picture, built as a perturbation around a trivial limit, can provide a highly accurate 
quantitative description of the system, when carried out to high orders. 
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TABLE I. Series coefficients for the integrated structure factor S(k) = m a n>m \ m cos(nfcd/2). Nonzero coefficients a. 



up to order m — 13 are listed. 



(n, m) 


0. n ,m 




(n, m) 






(n, m) 


a n m 




(n, m) 




1 

"3— 1 


( o, 0) 


1.000000000 




( 4,11) 


-1.591179666x10" 


3 


( 9, 7) 


-2.661559650x10 


-2 


(15, 9) 


-9.209089311x10" 




( !> 0) 


-1.000000000 




( 4,12) 


-1.346799079x10" 


3 


( 9, 8) 


-1.701711998x10" 


-2 


(15,10) 


-1.132484677x10" 


2 


( !. !) 


-2.500000000x10" 


l 


( 4,13) 


-1.135073475x10" 


3 


( 9, 9) 


-1.047345923x10 


-2 


(15,11) 


-1.152471081x10" 


2 


( !. 2) 


3.125000000x10" 


2 


( 5, 2) 


-9.375000000x10" 


2 


( 9,10) 


-6.243270719x10 


-3 


(15,12) 


-1.055798426x10" 


2 


( !> 3) 


2.864583333x10" 


2 


( 5, 3) 


-1.119791667x10" 


1 


( 9,11) 


-3.511373171x10 


-3 


(15,13) 


-9.126658607x10" 


3 


( !. 4 ) 


1.247829861x10" 


3 


( 5, 4) 


-5.577256944x10" 


2 


( 9,12) 


-1.832316909x10 


-3 


(16, 8) 


1.534223557x10" 


3 


( !. 5) 


4.747178819x10" 


4 


( 5, 5) 


-1.284450955x10" 


2 


( 9,13) 


-8.384814283x10 


-4 


(16, 9) 


4.630391871x10" 


3 


( !> 6) 


2.256346338x10" 


3 


( 5, 6) 


1.543539542x10" 


4 


(10, 5) 


1.538085938x10 


-2 


(16,10) 


7.630811251x10" 


3 


( !. 7) 


1.115188285x10" 


3 


( 5, 7) 


2.104722325x10" 


3 


(10, 6) 


2.759165823x10 


-2 


(16,11) 


9.315049401x10" 


3 


( !. 8) 


4.433909083x10" 


4 


( 5, 8) 


2.590777271x10" 


3 


(10, 7) 


2.762316480x10 


-2 


(16,12) 


9.630494889x10" 


3 


( !. 9) 


4.463708690x10" 


4 


( 5, 9) 


2.541851515x10" 


3 


(10, 8) 


2.143434446x10 


-2 


(16,13) 


9.064708700x10" 


3 


( 1,10) 


3.908585581x10" 


4 


( 5,10) 


2.142535125x10" 


3 


(10, 9) 


1.509559320x10 


-2 


(17, 8) 


-7.671117783x10" 


4 


( 1,H) 


2.721746108x10" 


4 


( 5,11) 


1.784799808x10" 


3 


(10,10) 


1.026773209x10 


-2 


(17, 9) 


-2.828093400x10" 


3 


( 1,12) 


2.078935750x10" 


4 


( 5,12) 


1.526601703x10" 


3 


(10,11) 


6.786844453x10 


-3 


(17,10) 


-5.469609444x10" 


3 


( 1,13) 


1.752731379x10" 


4 


( 5,13) 


1.307569342x10" 


3 


(10,12) 


4.357777095x10 


-3 


(17,11) 


-7.545021498x10" 


3 


( 2 , 1) 


5.000000000x10" 


1 


( 6, 3) 


7.812500000x10" 


2 


(10,13) 


2.725607152x10 


-3 


(17,12) 


-8.525536617x10" 


3 


( 2, 2) 


6.250000000x10" 


2 


( 6, 4) 


7.595486111x10" 


2 


(11, 5) 


-7.690429688x10 


-3 


(17,13) 


-8.551663397x10" 


3 


( 2, 3) 


-3.645833333x10" 


2 


( 6, 5) 


3.919813368x10" 


2 


(11, 6) 


-1.887210799x10 


-2 


(18, 9) 


7.244944572x10" 


4 


( 2, 4) 


-1.443142361x10" 


2 


( 6, 6) 


1.485306540x10" 


2 


(11, 7) 


-2.397200110x10 


-2 


(18,10) 


2.481638461x10" 


3 


( 2, 5) 


-3.861038773x10" 


3 


( 6, 7) 


5.066873213x10" 


3 


(11, 8) 


-2.190082988x10 


-2 


(18,11) 


4.609868075x10" 


3 


( 2, 6) 


-4.464373176x10" 


3 


( 6, 8) 


1.200659217x10" 


3 


(11, 9) 


-1.703925617x10 


-2 


(18,12) 


6.277698444x10" 


3 


( 2, 7) 


-3.574776355x10" 


3 


( 6, 9) 


-4.880865156x10" 


4 


(11,10) 


-1.243639951x10 


-2 


(18,13) 


7.142219966x10" 


3 


( 2, 8) 


-2.152545002x10" 


3 


( 6,10) 


-1.061951430x10" 


3 


(11,11) 


-8.818602901x10 


-3 


(19, 9) 


-3.622472286x10" 


4 


( 2, 9) 


-1.614262178x10" 


3 


( 6,11) 


-1.168436483x10" 


3 


(11,12) 


-6.103055681x10 


-3 


(19,10) 


-1.483667264x10" 


3 


( 2,10) 


-1.364268756x10" 


3 


( 6,12) 


-1.145143702x10" 


3 


(11,13) 


-4.131112380x10 


-3 


(19,11) 


-3.187079773x10" 


3 


( 2,11) 


-1.085128141x10" 


3 


( 6,13) 


-1.075193082x10" 


3 


(12, 6) 


7.049560547x10 


-3 


(19,12) 


-4.863723272x10" 


3 


( 2,12) 


-8.633290964x10" 


4 


( 7, 3) 


-3.906250000x10" 


2 


(12, 7) 


1.552691676x10 


-2 


(19,13) 


-6.030015201x10" 


3 


( 2,13) 


-7.178656848x10" 


4 


( 7, 4) 


-6.331380208x10" 


2 


(12, 8) 


1.896538987x10 


-2 


(20,10) 


3.441348672x10" 


4 


( 3, 1) 


-2.500000000x10" 


1 


( 7, 5) 


-4.956506800x10" 


2 


(12, 9) 


1.764401420x10 


-2 


(20,11) 


1.318864692x10" 


3 


( 3, 2) 


-1.875000000x10" 


1 


( 7, 6) 


-2.639009923x10" 


2 


(12,10) 


1.446797193x10 


-2 


(20,12) 


2.724516711x10" 


3 


( 3, 3) 


-2.343750000x10" 


2 


( 7, 7) 


-1.199321708x10" 


2 


(12,11) 


1.121713674x10 


-2 


(20,13) 


4.090782507x10" 


3 


( 3, 4) 


2.012803819x10" 


2 


( 7, 8) 


— 5.185196227X 10" 


3 


(12,12) 


8.435888979x10" 


-3 


(21,10) 


— 1.720674336X 10" 


4 


( 3, 5) 


1.035789207x10" 


2 


( 7, 9) 


-1.776089397x10" 


3 


(12,13) 


6.212068934x10 


-3 


(21,11) 


-7.750395765x10" 


4 


( 3, 6) 


4.805176346x10" 


3 


( 7,10) 


-8.192333266x10" 


5 


(13, 6) 


-3.524780273x10 


-3 


(21,12) 


-1.829266700x10" 


3 


( 3, 7) 


4.284713101x10" 


3 


( 7,11) 


6.216205777x10" 


4 


(13, 7) 


-1.010268412x10 


-2 


(21,13) 


-3.055890338x10" 


3 


( 3, 8) 


3.284367466x10" 


3 


( 7,12) 


8.737107987x10" 


4 


(13, 8) 


-1.511986766x10 


-2 


(22,11) 


1.642461866x10" 


4 


( 3, 9) 


2.260859277x10" 


3 


( 7,13) 


9.485913165x10" 


4 


(13, 9) 


-1.626035314x10 


-2 


(22,12) 


6.962934102x10" 


4 


( 3,10) 


1.762829643x10" 


3 


( 8, 4) 


3.417968750x10" 


2 


(13,10) 


-1.465683078x10 


-2 


(22,13) 


1.582696489x10" 


3 


( 3,11) 


1.451718673x10" 


3 


( 8, 5) 


4.730902778x10" 


2 


(13,11) 


-1.210970065x10 


-2 


(23,11) 


-8.212309331x10" 


5 


( 3,12) 


1.173179929x10" 


3 


( 8, 6) 


3.645070394x10" 


2 


(13,12) 


-9.584599362x10 


-3 


(23,12) 


-4.034250937x10" 


4 


( 3,13) 


9.595019662x10" 


4 


( 8, 7) 


2.204519358x10" 


2 


(13,13) 


-7.395426977x10 


-3 


(23,13) 


-1.037148415x10" 


3 


( 4, 2) 


1.875000000x10" 


1 


( 8, 8) 


1.249375691x10" 


2 


(14, 7) 


3.273010254x10 


-3 


(24,12) 


7.870129775x10" 




( 4, 3) 


1.041666667x10" 


1 


( 8, 9) 


6.866897220x10" 


3 


(14, 8) 


8.544238997x10" 


-3 


(24,13) 


3.656570992x10" 


4 


( 4, 4) 


1.909722222x10" 


2 


( 8,10) 


3.496121257x10" 


3 


(14, 9) 


1.227229300x10 


-2 


(25,12) 


-3.935064888x10" 


5 


( 4, 5) 


-3.906250000x10" 


3 


( 8,11) 


1.579388581x10" 


3 


(14,10) 


1.322238413x10 


-2 


(25,13) 


-2.093578942x10" 


4 


( 4, 6) 


-3.848888256x10" 


3 


( 8,12) 


5.538137269x10" 


4 


(14,11) 


1.227668397x10 


-2 


(26,13) 


3.783716238x10" 


5 


( 4, 7) 


-3.145002043x10" 


3 


( 8,13) 


4.738508965x10" 


6 


(14,12) 


1.055795986x10 


-2 


(27,13) 


-1.891858119x10" 


5 


( 4, 8) 


-2.985769420x10" 


3 


( 9, 4) 


-1.708984375x10" 


2 


(14,13) 


8.705553180x10 


-3 








( 4, 9) 


-2.431829046x10" 


3 


( 9, 5) 


-3.485333478x10" 


2 


(15, 7) 


-1.636505127x10 


-3 








( 4,10) 


-1.912181882x10" 


3 


( 9, 6) 


-3.606061582x10" 


2 


(15, 8) 


-5.362708709x10 


-3 
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TABLE II. Series coefficients for the exclusive structure factors of the 1-particle triplet state Si p (k) Si p (k) = m a. 



A m cos(nfcd/2). Nonzero coefficients a n , m up to order m = 13 are listed. 



(n, m) 






(n, m) 






(n, m) 


an, ?7i 




(n, m) 






( 0, 0) 


1.000000000 




( 4. 2) 


1.875000000x10" 


1 


( 8,11) 


-6.011492935x10" 


3 


(14,13) 


8.674883160x10" 




( 0, 2) 


-3.125000000x10" 


1 


( 4, 3) 


1.458333333x10" 


1 


( 8,12) 


9.711804995x10" 


3 


(15, 7) 


-1.636505127x10" 




( 0, 3) 


-9.375000000x10" 


2 


( 4, 4) 


6.835937500x10" 


3 


( 8,13) 


4.710318350x10" 


3 


(15, 8) 


-5.358020138x10" 




( 0, 4) 


1.627604167x10" 


2 


( 4, 5) 


-5.485930266x10" 


2 


( 9, 4) 


-1.708984375x10" 


2 


(15, 9) 


-9.241334017x10" 


3 


( 0, 5) 


-8.257378472x10" 


2 


( 4, 6) 


9.334422924x10" 


3 


( 9, 5) 


-3.543203848x10" 


2 


(15,10) 


-1.131045012x10" 


2 


( 0, 6) 


-4.077148438x10" 


2 


( 4, 7) 


7.197318057x10" 


3 


( 9, 6) 


-3.456574899x10" 


2 


(15,11) 


-1.142877596x10" 


2 


( 0, 7) 


1.841892038x10" 


2 


( 4, 8) 


-5.832889946x10" 


2 


( 9, 7) 


-2.518777141x10" 


2 


(15,12) 


-1.046099406x10" 


2 


( 0, 8) 


-4.383319893x10" 


2 


( 4, 9) 


-1.281724218x10" 


2 


( 9, 8) 


-1.987998380x10" 


2 


(15,13) 


-9.327327593x10" 


3 


( 0, 9) 


-3.392954635x10" 


2 


( 4,10) 


2.152346878x10" 


2 


( 9, 9) 


-1.377372682x10" 


2 


(16, 8) 


1.534223557x10" 


3 


( 0,10) 


1.671238124x10" 


2 


( 4,11) 


-4.067856177x10" 


2 


( 9,10) 


-1.514186704x10" 


3 


(16, 9) 


4.632624525x10" 


3 


( 0,11) 


-2.402930407x10" 


2 


( 4,12) 


-2.401466656x10" 


2 


( 9,11) 


-2.640524172x10" 


4 


(16,10) 


7.616553597x10" 


3 


( 0,12) 


-3.094915187x10" 


2 


( 4,13) 


2.523955463x10" 


2 


( 9,12) 


-8.667418864x10" 


3 


(16,11) 


9.310862860x10" 


3 


( 0,13) 


1.317777693x10" 


2 


( 5, 2) 


-9.375000000x10" 


2 


( 9,13) 


-3.106212436x10" 


4 


(16,12) 


9.677973448x10" 


3 


( !. 0) 


-1.000000000 




( 5, 3) 


-1.223958333x10" 


1 


(10, 5) 


1.538085938x10" 


2 


(16,13) 


9.174113955x10" 


3 


( 1> 1) 


-2.500000000x10" 


1 


( 5, 4) 


-4.031032986x10" 


2 


(10, 6) 


2.725408107x10" 


2 


(17, 8) 


-7.671117783x10" 


4 


( 1> 2 ) 


5.000000000x10" 


1 


( 5, 5) 


-1.652470341x10" 


3 


(10, 7) 


2.832565779x10" 


2 


(17, 9) 


-2.828986461x10" 


3 


( !. 3) 


1.432291667x10" 


1 


( 5, 6) 


-2.680371131x10" 


2 


(10, 8) 


2.311140563x10" 


2 


(17,10) 


-5.461641059x10" 


3 


( !. 4 ) 


-3.797743056x10" 


2 


( 5, 7) 


-2.120727570x10" 


3 


(10, 9) 


1.449731923x10" 


2 


(17,11) 


-7.554290412x10" 


3 


( 1> 5) 


1.475604022x10" 


1 


( 5, 8) 


3.677335722x10" 


2 


(10,10) 


6.246098908x10" 


3 


(17,12) 


-8.549734505x10" 


3 


( !. 6) 


6.164414206x10" 


2 


( 5, 9) 


-4.177481543x10" 


3 


(10,11) 


6.405727927x10" 


3 


(17,13) 


-8.565778023x10" 


3 


( 1. 7) 


-4.571135603x10" 


2 


( 5,10) 


-1.510927532x10" 


2 


(10,12) 


8.071629661x10" 


3 


(18, 9) 


7.244944572x10" 


4 


( 1. 8) 


7.541900211x10" 


2 


( 5.H) 


3.391969899x10" 


2 


(10,13) 


-1.769051814x10" 


4 


(18,10) 


2.481229141x10" 


3 


( 1. 9) 


5.347557773x10" 


2 


( 5,12) 


1.090426758x10" 


2 


(11, 5) 


-7.690429688x10" 


3 


(18,11) 


4.613304272x10" 


3 


( 1,10) 


-4.190407427x10" 


2 


( 5,13) 


-2.101789738x10" 


2 


(11, 6) 


-1.875154472x10" 


2 


(18,12) 


6.275948659x10" 


3 


( l.n) 


3.910767511x10" 


2 


( 6, 3) 


7.812500000x10" 


2 


(11, 7) 


-2.441788701x10" 


2 


(18,13) 


7.129091829x10" 




( 1.12) 


5.118486369x10" 


2 


( 6, 4) 


6.727430556x10" 


2 


(11, 8) 


-2.218640450x10" 


2 


(19, 9) 


-3.622472286x10" 


4 


( 1.13) 


-3.356079789x10" 


2 


( 6, 5) 


4.458279080x10" 


2 


(11, 9) 


-1.605467746x10" 


2 


(19,10) 


-1.483499815x10" 




( 2 . 1) 


5.000000000x10" 


1 


( 6, 6) 


3.504162070x10" 


2 


(11,10) 


-1.111541180x10" 


2 


(19,11) 


-3.188966289x10" 


3 


( 2, 2) 


-1.250000000x10" 


1 


( 6, 7) 


2.109803094x10" 


3 


(11,11) 


-1.034165700x10" 


2 


(19,12) 


-4.859982468x10" 


3 


( 2, 3) 


-2.604166667x10" 


2 


( 6, 8) 


-1.579822777x10" 


2 


(11,12) 


-8.030332288x10" 


3 


(19,13) 


-6.025154806x10" 


3 


( 2, 4) 


1.779513889x10" 


2 


( 6, 9) 


1.297988287x10" 


2 


(11,13) 


-1.448277762x10" 


3 


(20,10) 


3.441348672x10" 


4 


( 2, 5) 


-1.206416377x10" 


1 


( 6,10) 


6.764275332x10" 


3 


(12, 6) 


7.049560547x10" 


3 


(20,11) 


1.318939114x10" 


3 


( 2, 6) 


-2.833321654x10" 


2 


( 6,11) 


-2.652577291x10" 


2 


(12, 7) 


1.559121717x10" 


2 


(20,12) 


2.723720941x10" 


3 


( 2, 7) 


4.500562644x10" 


2 


( 6,12) 


-1.873115433x10" 


3 


(12, 8) 


1.875765450x10" 


2 


(20,13) 


4.091885849x10" 


3 


( 2, 8) 


-6.561594063x10" 


2 


( 6,13) 


1.362791590x10" 


2 


(12, 9) 


1.728643232x10" 


2 


(21,10) 


-1.720674336x10" 


4 


( 2, 9) 


-3.427604879x10" 


2 


( 7, 3) 


-3.906250000x10" 


2 


(12,10) 


1.481153271x10" 


2 


(21,11) 


-7.750705856x10" 


4 


( 2,10) 


4.545736434x10" 


2 


( 7, 4) 


-6.070963542x10" 


2 


(12,11) 


1.267199966x10" 


2 


(21,12) 


-1.828835480x10" 


3 


( 2 .H) 


-3.409665942x10" 


2 


( 7, 5) 


-5.429642289x10" 


2 


(12,12) 


8.696286848x10" 


3 


(21,13) 


-3.057141912x10" 


3 


( 2,12) 


-3.879019409x10" 


2 


( 7, 6) 


-3.120500070x10" 


2 


(12,13) 


4.298281913x10" 


3 


(22,11) 


1.642461866x10" 


4 


( 2,13) 


3.772049749x10" 


2 


( 7, 7) 


-3.188448384x10" 


3 


(13, 6) 


-3.524780273x10" 


3 


(22,12) 


6.962799729x10" 


4 


( 3, 1) 


-2.500000000x10" 


1 


( 7, 8) 


1.279329754x10" 


3 


(13, 7) 


-1.012679677x10" 


2 


(22,13) 


1.582874937x10" 


3 


( 3, 2) 


-1.562500000x10" 


1 


( 7, 9) 


-1.467444993x10" 


2 


(13, 8) 


-1.499607686x10" 


2 


(23,11) 


-8.212309331x10" 




( 3, 3) 


-8.593750000x10" 


2 


( 7,10) 


-9.349685042x10" 


4 


(13, 9) 


-1.624052105x10" 


2 


(23,12) 


-4.034194087x10" 


4 


( 3, 4) 


1.372612847x10" 


2 


( 7,11) 


1.714891750x10" 


2 


(13,10) 


-1.498198003x10" 


2 


(23,13) 


-1.037244158x10" 


3 


( 3, 5) 


1.005768953x10" 


1 


( 7,12) 


-5.294607061x10" 


3 


(13,11) 


-1.252169240x10" 


2 


(24,12) 


7.870129775x10" 




( 3, 6) 


9.286244710x10" 


3 


( 7,13) 


-8.581846951x10" 


3 


(13,12) 


-9.049303257x10" 


3 


(24,13) 


3.656595110x10" 


4 


( 3, 7) 


-2.312368150x10" 


2 


( 8. 4) 


3.417968750x10" 


2 


(13,13) 


-6.457204366x10" 


3 


(25,12) 


-3.935064888x10" 


5 


( 3, 8) 


6.813065057x10" 


2 


( 8, 5) 


4.904513889x10" 


2 


(14, 7) 


3.273010254x10" 


3 


(25,13) 


-2.093589278x10" 


4 


( 3, 9) 


2.664720824x10" 


2 


( 8, 6) 


3.434541490x10" 


2 


(14, 8) 


8.532182670x10" 


3 


(26,13) 


3.783716238x10" 


5 


( 3,10) 


-3.218290488x10" 


2 


( 8, 7) 


1.559162061x10" 


2 


(14, 9) 


1.232853596x10" 


2 


(27,13) 


-1.891858119x10" 


5 


( 3,11) 


4.070336366x10" 


2 


( 8, 8) 


1.322605786x10" 


2 


(14,10) 


1.328042130x10" 


2 








( 3,12) 


3.466875907x10" 


2 


( 8, 9) 


1.580418652x10" 


2 


(14,11) 


1.213368398x10" 


2 








( 3,13) 


-3.003621683x10" 


2 


( 8,10) 


9.329997203x10" 


4 


(14,12) 


1.012086983x10" 


2 
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TABLE III. Series coefficients for the integrated structure factor S(k), the exclusive 1-particle structure factor S\ p (k), and 
the total 2-particle structure factor S2 P (k) at kd = it, 2tv, together with the structure factors (Sti(&), ST 2 {k), and S2 P c(k)) for 
2-particle bound states T\ and T2, and 2-particle continuum at kd = it, the auto correlation function of 1-particle ($i p ) and 
2-particle states (<&2 P ), the average relative weight for 2-particle states (W2 P ), and R, Ri p , and 7?2 P . Series coefficients of \ m 
up to order m = 13 are listed. 



m 


SO) 








S2pW 




St, M 






1 





1.000000000 




1.000000000 




0.000000000 




0.000000000 




0.000000000 




1 


-5.000000000x10" 




-5.000000000x10" 




0.000000000 




0.000000000 




0.000000000 




2 


1.250000000x10" 




0.000000000 




1.250000000x10" 




1.250000000x10" 




0.000000000 




3 


6.250000000x10" 




0.000000000 




6.250000000x10" 




6.250000000x10" 




0.000000000 




4 


-8.246527778x10" 




-2.777777778x10" 




1.584201389x10" 




4.557291667x10" 




2.604166667x10" 




5 


-7.315176505x10" 




-2.770996094x10" 




1.725260417x10" 




-2.365451389x10" 




5.642361111x10" 




6 


l.ur lUzo/ozXIU 


3 


— z.4UU4o ( IzoX 10 


2 


z.ooi iz4yoo x io 


2 


— 6.690o094oz X 10 


2 


4.o4ol9z9yo X 10 


3 


7 


2.038836381x10" 


3 


-2.191502135x10" 


2 


2.248861054x10" 


2 


-1.295114399x10" 


1 


-3.105917095x10" 


3 


8 


9.809032520 x 10" 


4 


— 1.887358236X 10" 


2 


1.753652230x10" 


2 


-1.919116349x10" 


1 


-1.721538889x10" 


2 


9 


7. 194422822 x 10" 


4 


— 1.527772889X 10" 


2 


1.378672204x10" 


2 


-2.252182500x10" 


1 


-3.106542706x10" 


2 


If) 


4.813229257 x 10" 


4 


— 1.22883 1810 x 10" 


2 


1.083241128x10" 


2 


-1.979558323x10" 


1 


-3.789373385x10" 


2 


11 


2.551816973 x 10" 


4 


— 1.011208718 x 10" 


2 


9.645537754x10" 


3 


-7.313449239x10" 


2 


-3.170664884x10" 


2 


12 


1. 953605167 x 10" 


4 


— 8. 576749497 x 10" 


3 


7.501091661x10" 


3 


1.758607061x10" 


1 


-6.069523293x10" 


3 


13 


2 090270912x10" 


4 


— 7 538604157x10" 


3 


7.444260493x10" 


3 


5.577962015x10" 


1 


4.352098881x10" 


2 


14 


















1.144814294x10" 


1 


m 


5 2 pc(7T) 




o {2tt ) 




Si p (27r) 




S 2p (27r) 






1 





u.uuuuuuuuu 




z.UUUUUUUU U 




2.000000000 




0.000000000 




1.000000000 




1 


n nnnnnnnnn 

U . U \-J U U U U U U U 




i nnnnnnnnn 




1.000000000 




0.000000000 




0.000000000 




2 


n nnnnnnnnn 
u . u u u u u u u u u 




^ nnnnnnnnn v i n ~ 


1 


-5.000000000x10" 


1 


1.000000000 




-3.125000000x10" 


1 


3 


n nnnnnnnnn 
u . u u u u u u u u u 




9 Q1fifif!fififi7Yin" 
Z.yiODDDDO I X ±U 


1 


2.083333333x10" 


1 


8.333333333x10" 


2 


-9.375000000x10" 


2 


4 


O.DoUOOOOOD X _LU 


3 


9 9QfinnfiQzLA v i n ~ 


1 


2.847222222x10" 


1 


-9.027777778x10" 


2 


1.627604167x10" 


2 


5 


3.526475694x10" 


2 


1.882414641x10" 


1 


-2.981318721x10" 


1 


4.924768519x10" 


1 


-8.257378472x10" 


2 


6 


8.583615150x10" 


2 


1.552634534x10" 


1 


8.784079846x10" 


2 


1.213981723x10" 


1 


-4.077148438x10" 


2 


7 


1.551059676x10" 


1 


1.336307604x10" 


1 


2.710263476x10" 


1 


-1.130861384x10" 


1 


1.841892038x10" 


2 


8 


2.266635461x10" 


1 


1.180685972x10" 


1 


-2.368294851x10" 


1 


3.319204641x10" 


1 


-4.383319893x10" 


2 


9 


2.700703991x10" 


1 


1.053990124x10" 


1 


-5.538722908x10" 


3 


2.095799705x10" 


1 


-3.392954635x10" 


2 


10 


2.466819774x10" 


1 


9.514478385x10" 


2 


2.723409199x10" 


1 


3.264975567x10" 


2 


1.671238124x10" 


2 


11 


1.144866790x10" 


1 


8.684667562x10" 


2 


-1.694460542x10" 


1 


4.627752009x10" 


1 


-2.402930407x10" 


2 


12 


-1.622900912x10" 


1 


7.991574506x10" 


2 


-7.914782458x10" 


2 


4.466072907x10" 


1 


-3.094915187x10" 


2 


13 


-5.938729299x10" 


1 


7.400747491x10" 


2 


2.593075729x10" 


1 


7.243209793x10" 


1 


1.317777693x10" 


2 


m 










R = \im k ^ S/k 2 


Ri v = lim fc ^o Sip/k 2 


R2 V = lim fc ^o S-2-p/k 2 \ 





0.000000000 




1.000000000 




1.250000000x10" 


1 


1.250000000x10" 


1 


0.000000000 




1 


0.000000000 




0.000000000 




6.250000000x10" 


2 


6.250000000x10" 


2 


0.000000000 




2 


3.125000000x10" 


1 


-1.875000000x10" 


1 


9.375000000x10" 


2 


9.375000000x10" 


2 


0.000000000 




3 


9.375000000x10" 


2 


-4.687500000x10" 


2 


7.031250000x10" 


2 


7.031250000x10" 


2 


0.000000000 




4 


-2.636718750x10" 


2 


-3.103298611x10" 


2 


6.610785590x10" 


2 


6.141493056x10" 


2 


4.448784722x10" 


3 


5 


8.148871528x10" 


2 


-6.538447627x10" 


2 


6.383938260x10" 


2 


6.599765354x10" 


2 


-2.305772569x10" 


3 


6 


5.053823966x10" 


2 


-1.166864089x10" 


2 


6.018220642x10" 


2 


6.027503661x10" 


2 


-1.646442178x10" 


4 


7 


-1.360771980x10" 


2 


2.293814357x10" 


3 


5.737893768x10" 


2 


5.616068546x10" 


2 


9.660760071x10" 


4 


8 


3.967389006x10" 


2 


-3.920166318x10" 


2 


5.536595911x10" 


2 


5.331091464x10" 


2 


1.625177465x10" 


3 


9 


5.235992915x10" 


2 


-1.075891724x10" 


2 


5.357448621x10" 


2 


5.601825866x10" 


2 


-3.020253684x10" 


3 


10 


1.907595256x10" 


2 


7.139608316x10" 


3 


5.194887667x10" 


2 


5.101011394x10" 


2 


-3.106958042x10" 


5 


11 


6.149786123x10" 


2 


-2.797421962x10" 


2 


5.053076987x10" 


2 


4.803748655x10" 


2 


6.937664269x10" 


4 


12 


7.568512290x10" 


2 


-1.386098936x10" 


2 


4.927953902x10" 


2 


5.120296763x10" 


2 


-5.289988311x10" 


3 


13 


1.267526913x10" 


1 


9.053856936x10" 


3 


4.815020988x10" 


2 


4.831429881x10" 


2 


-6.876775997x10" 


3 
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TABLE IV. [n/m] Dlog Pade approximants to the series for dS 2/3 /d\ and 8S 4/5 /dX at kd = 2tt, for $i p , and for d$ 2p /d\. 
The position of the singularity (pole), the critical index (U.B.R.) from unbiased approximants, and the critical index from 
biased approximants (B.R.) are given. An asterisk denotes a defective approximant. 



n [(n-2)/n\ [(n - l)/n] [n/n] [(n + l)/n] [(n + 2)/n] 

polo (N.B.R. , B.R.) polo (N.B.R., B.R.) pole (N.B.R., B.R.) polo (N.B.R., B.R.) polo (N.B.R., B.R.) 

dS(2TT) 2 ' 3 /d\ 

n= 1 0.6244(-0.2599, -0.7844) 1.0859(-1.3670, -0.9980) 1.0837(-1.3560, -0.2298) 

n= 2 * ( * -1.0983) * ( * , * ) 1.0837(-1.3561, -0.9488) * ( * , -0.9530) 0.9694(-0.7956, -0.9621) 

n= 3 * ( * -1.1246) 0.9952(-0.9312, -0.9534) 0.9895(-0.9048, -0.9449) 0.9933(-0.9257, -0.9802) 0.9974(-0.9523, -0.9747) 

n= 4 0.9899(-0.9074, -0.9629) 0.9918(-0.9168, -0.9876) * ( * . -0.9738) 0.9981(-0.9582, -0.9765) 0.9980(-0.9576, -0.9300) 

n= 5 * ( * -0.9764) 0.9970(-0.9483, -0.9773) 0.9980(-0.9575, -0.9821) 0.9981(-0.9581, -0.9832) 

n= 6 0.9981(-0.9586, -0.9758) 0.9993(-0.9724, -0.9832) 

dS(2n) 4/5 /d\ 

n= 1 0.6406(-0.2873, -0.8325) 1.0777(-1.3676,-1.0271) 1.0833(-1.3962, * ) 

n= 2 * ( * , -1.0775) 0.8684(-0.6358,-1.2702) 1.0833(-1.3967, -0.9784) 1.0779(-1.3683,-0.9813) 0.9710(-0.8264, -0.9895) 

n= 3 * ( * -1.1123) 0.9969(-0. 9668,-0. 9814) 0.9907(-0.9375, -0.9768) 0.9941(-0. 9565,-1. 0057) 0.9979(-0.9823, -1.0000) 

n= 4 0.9913(-0.9406, -0.9898) 0.9930(-0. 9494,-1. 0134) * ( * . -0.9990) 0.9986(-0. 9875,-1. 0015) 0.9985(-0.9865, -0.9861) 

n= 5 * ( * -1.0016) 0.9975(-0. 9782,-1. 0021) 0.9984(-0.9862, -1.0058) 0.9985(-0. 9872,-1. 0065) 

n= 6 0.9985(-0.9871, -1.0014) 0.9995(-0. 9994.-1. 0065) 

n= 1 * ( * 0.0118) 0.2328(0.0004, 0.2419) * ( * , 0.3983) 

n= 2 * ( * 0.2984) * ( * . 0.3254) 0.6947(0.0613, 0.3328) * ( * > 0.3252) * ( * . 0.3231) 

n= 3 0.9773(0.3025, 0.3354) 1.0250(0.3680. 0.3289) 1.0065(0.3359, 0.3224) 0.9901(0.3043, 0.3244) 0.9935(0.3117, 0.3327) 

n= 4 1.0111(0.3449, 0.3250) 0.9266(0.1528. 0.3257) 0.9929(0.3103, 0.3269) * ( * , 0.3311) 0.9953(0.3162, 0.3323) 

n= 5 0.9965(0.3188, 0.3215) 0.9954(0.3164, 0.3218) 0.9948(0.3148, 0.3187) 0.9952(0.3158, 0.3568) 

n= 6 * ( * , 0.3215) 0.9951(0.3154, 0.3225) 

9*2 P /9A 

n= 1 -0.2374( 0.0304, -0.1924) * ( * , -1.1986) * ( * * ) 

n= 2 0.7175(-0.1129, -0.2521) 0.2556(-0.0140, -0.1334) ( * , -0.7531) * ( * , -0.7601) * ( * * ) 

n= 3 0.4363(-0.0604, 0.0112) 0.9908(-0.7323, -0.7603) 0.7409(-0.2107, -0.7529) 0.6647(-0.1138, -0.4809) 0.5793(-0.0453, -0.6359) 

n= 4 0.8362(-0.3880, -1.0869) 0.6399(-0.0869. -0.6158) 0.3412(-0.0003, -0.6030) 0.6184(-0.0724, -0.5041) 0.6071(-0.0630, -0.1906) 

n= 5 0.4558(-0.0054, -0.6038) 0.5885(-0.0485, -0.6162) 0.6023(-0.0588, * ) 
n= 6 0.6035(-0.0599, -0.4355) 



2G 



